Candidate Genes Research Articles

Multiomics-assisted characterization of rice-Yellow Stem Borer interaction provides genomic and mechanistic insights into stem borer resistance in rice.

By deploying a multi-omics approach, we unraveled the mechanisms that might help rice to combat Yellow Stem Borer infestation, thus providing insights and scope for developing YSB resistant rice varieties. Yellow Stem Borer (YSB), Scirpophaga incertulas (Walker) (Lepidoptera: Crambidae), is a major pest of rice, that can lead to 20-60% loss in rice production. Effective management of YSB infestation is challenged by the non-availability of adequate sources of resistance and poor understanding of resistance mechanisms, thus necessitating studies for generating resources to breed YSB resistant rice and to understand rice-YSB interaction. In this study, by using bulk-segregant analysis in combination with next-generation sequencing, Quantitative Trait Loci (QTL) intervals in five rice chromosomes were mapped that could be associated with YSB resistance at the vegetative phase in a resistant rice line named SM92. Further, multiple SNP markers that showed significant association with YSB resistance in rice chromosomes 1, 5, 10, and 12 were developed. RNA-sequencing of the susceptible and resistant lines revealedseveral genes present in the candidate QTL intervals to be differentially regulated upon YSB infestation. Comparative transcriptome analysis revealed a putative candidate gene that was predicted to encode an alpha-amylase inhibitor. Analysis of the transcriptome and metabolite profiles further revealed a possible link between phenylpropanoid metabolism and YSB resistance. Taken together, our study provides deeper insights into rice-YSB interaction and enhances the understanding of YSB resistance mechanism. Importantly, a promising breeding line and markers for YSB resistance have been developed that can potentially aid in marker-assisted breeding of YSB resistance among elite rice cultivars.

Intragenic MFSD8 duplication and histopathological findings in a rabbit with neuronal ceroid lipofuscinosis.

Neuronal ceroid lipofuscinoses (NCL) are among the most prevalent neurodegenerative disorders of early life in humans. Disease-causing variants have been described for 13 different NCL genes. In this study, a refined pathological characterization of a female rabbit with progressive neurological signs reminiscent of NCL was performed. Cytoplasmic pigment present in neurons was weakly positive with Sudan black B and autofluorescent. Immunohistology revealed astrogliosis, microgliosis and axonal degeneration. During the subsequent genetic investigation, the genome of the affected rabbit was sequenced and examined for private variants in NCL candidate genes. The analysis revealed a homozygous ~10.7 kb genomic duplication on chromosome 15 comprising parts of the MFSD8 gene, NC_013683.1:g.103,727,963_103,738,667dup. The duplication harbors two internal protein coding exons and is predicted to introduce a premature stop codon into the transcript, truncating ~50% of the wild-type MFSD8 open reading frame encoding the major facilitator superfamily domain containing protein 8, XP_002717309.2:p.(Glu235Leufs*23). Biallelic loss-of-function variants in MFSD8 have been described to cause NCL7 in human patients, dogs and a single cat. The available clinical and pathological data, together with current knowledge about MFSD8 variants and their functional impact in other species, point to the MFSD8 duplication as a likely causative defect for the observed phenotype in the affected rabbit.

Genetic causes and genomic consequences of breakdown of distyly in Linum trigynum.

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility (SI). Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male SI and/or anther height could underlie homostyly and self-compatibility (SC) in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.

Novel candidate genes for environmental stresses response in Synechocystis sp. PCC 6803 revealed by machine learning algorithms.

Cyanobacteria have developed acclimation strategies to adapt to harsh environments, making them a model organism. Understanding the molecular mechanisms of tolerance to abiotic stresses can help elucidate how cells change their gene expression patterns in response to stress. Recent advances in sequencing techniques and bioinformatics analysis methods have led to the discovery of many genes involved in stress response in organisms. The Synechocystis sp. PCC 6803 is a suitable microorganism for studying transcriptome response under environmental stress. Therefore, for the first time, we employed two effective feature selection techniques namely and support vector machine recursive feature elimination (SVM-RFE) and LASSO (Least Absolute Shrinkage Selector Operator) to pinpoint the crucial genes responsive to environmental stresses in Synechocystis sp. PCC 6803. We applied these algorithms of machine learning to analyze the transcriptomic data of Synechocystis sp. PCC 6803 under distinct conditions, encompassing light, salt and iron stress conditions. Seven candidate genes namely sll1862, slr0650, sll0760, slr0091, ssl3044, slr1285, and slr1687 were selected by both LASSO and SVM-RFE algorithms. RNA-seq analysis was performed to validate the efficiency of our feature selection approach in selecting the most important genes. The RNA-seq analysis revealed significantly high expression for five genes namely sll1862, slr1687, ssl3044, slr1285, and slr0650 under ion stress condition. Among these five genes, ssl3044 and slr0650 could be introduced as new potential candidate genes for further confirmatory genetic studies, to determine their roles in their response to abiotic stresses.

Identification and map-based cloning of an EMS-induced mutation in wheat gene TaSP1 related to spike architecture.

A candidate gene TaSP1 related to spike shape was cloned, and the gene-specific marker was developed to efficiently track the superior haplotype in common wheat. Spike shape, an important factor that affects wheat grain yield, is mainly defined by spike length (SPL), spikelet number (SPN), and compactness. Zhoumai32 mutant 1160 (ZM1160), a mutant obtained from ethyl methane sulfonate (EMS) treatment of hexaploid wheat variety Zhoumai32, was used to identify and clone the candidate gene that conditioned the spike shape. Genetic analysis of an F2 population derived from a cross of ZM1160 and Bainong207 suggested that the compact spike shape in ZM1160 was controlled by a single recessive gene, and therefore, the mutated gene was designated as Tasp1. With polymorphic markers identified through bulked segregant analysis (BSA), the gene was mapped to a 2.65-cM interval flanked by markers YZU0852 and MIS46239 on chromosome 7D, corresponding to a 0.42-Mb physical interval of Chinese spring (CS) reference sequences (RefSeq v1.0). To fine map TaSP1, 15 and seven recombinants were, respectively, screened from 1599 and 1903 F3 plants derived from the heterozygous F2 plants. Finally, TaSP1 was delimited to a 21.9Kb (4,870,562 to 4,892,493bp) Xmis48123-Xmis48104 interval. Only one high-confidence gene TraesCS7D02G010200 was annotated in this region, which encodes an unknown protein with a putative vWA domain. Quantitative reverse transcription PCR (qRT-PCR) analysis showed that TraesCS7D02G010200 was mainly expressed in the spike. Haplotype analysis of 655 wheat cultivars using the candidate gene-specific marker Xg010200p2 identified a superior haplotype TaSP1b with longer spike and more spikelet number. TaSP1 is beneficial to the improvement in wheat spike shape.

Transcriptomic responses to drought stress in the Patagonian southern beech Nothofagus alpina

BackgroundDeciphering the genetic architecture of drought tolerance could allow the candidate genes identification responding to water stress. In the Andean Patagonian forest, the genus Nothofagus represents an ecologically relevant species to be included in different genomic studies. These studies are scarce in South American ecosystems however represent an important source of genomic data in order to interpret future climate-change environment scenarios of these emblematic forests. Here, we achieved the assemblage of the transcriptome of N. alpina while searching for key genes of activated or suppressed metabolic pathways in response to drought stress.ResultsDe novo transcriptome assembly resulted in 104,030 transcripts. Following confirmation of drought conditions, based on reduction of leaf water potential and stomatal conductance, a differential gene expression analysis resulted in 2720 significantly expressed genes (1601 up-regulated and 1119 down-regulated). Enrichment analysis (over-representation analysis and gene set enrichment analysis) resulted in more than one hundred stress-responsive term ontologies (i.e. biological processes) and pathways. Terms such as response to abscisic acid and pathways such as plant hormone signal transduction or starch and sucrose metabolism were over-represented. Protein–protein interaction assessment resulted in networks with significantly expressed top common hub gene clusters (e.g. plant-type cell wall biogenesis among down-regulated or ABA-signalling among up-regulated). These networks evidenced important regulators at gene expression such as transcriptional factors.ConclusionsResponses of N. alpina seedlings to drought stress were evidenced by the activation of several genes linked to GO biological processes and KEGG pathways, which were mainly based on over-expression of specific protein kinases, phosphatases, synthases and transcription factors. This suggests an up-regulation of signalling pathways, triggered through plant hormones such as abscisic acid or auxin, which could counteract the osmotic stress created as a probable immediate response to drought. On the other hand, groups of carbon fixation genes related to the galactose metabolism, photosynthesis, secondary wall biogenesis, and fatty acid biosynthesis degradation were down-regulated under drought. Overall, our results provide new genomic data for understanding how non-model long-lived trees of Patagonian forests would acclimate to environmental changes.

Identification of a candidate gene for the I locus determining the dominant white bulb color in onion (Allium cepa L.).

Through a map-based cloning approach, a gene coding for an R2R3-MYB transcription factor was identified as a causal gene for the I locus controlling the dominant white bulb color in onion. White bulb colors in onion (Allium cepa L.) are determined by either the C or I loci. The causal gene for the C locus was previously isolated, but the gene responsible for the I locus has not been identified yet. To identify candidate genes for the I locus, an approximately 7-Mb genomic DNA region harboring the I locus was obtained from onion and bunching onion (A. fistulosum) whole genome sequences using two tightly linked molecular markers. Within this interval, the AcMYB1 gene, known as a positive regulator of anthocyanin production, was identified. No polymorphic sequences were found between white and red AcMYB1 alleles in the 4,860-bp full-length genomic DNA sequences. However, a 4,838-bp LTR-retrotransposon was identified in the white allele, in the 79-bp upstream coding region from the stop codon. The insertion of this LTR-retrotransposon created a premature stop codon, resulting in the replacement of 26 amino acids with seven different residues. A molecular marker was developed based on the insertion of this LTR-retrotransposon to genotype the I locus. A perfect linkage between bulb color phenotypes and marker genotypes was observed among 5,303 individuals of segregating populations. The transcription of AcMYB1 appeared to be normal in both red and white onions, but the transcription of CHS-A, which encodes chalcone synthase and is involved in the first step of the anthocyanin biosynthesis pathway, was inactivated in the white onions. Taken together, an aberrant AcMYB1 protein produced from the mutant allele might be responsible for the dominant white bulb color in onions.

Fine-mapping and validation of the major quantitative trait locus QFlANG-4B for flag leaf angle in wheat.

This study precisely mapped and validated a quantitative trait locus (QTL) located on chromosome 4B for flag leaf angle in wheat. Flag leaf angle (FLANG) is closely related to crop architecture and yield. We previously identified the quantitative trait locus (QTL) QFLANG-4B for FLANG on chromosome 4B, located within a 14-cM interval flanked by the markers Xbarc20 and Xzyh357, using a mapping population of recombinant inbred lines (RILs) derived from a cross between Nongda3331 (ND3331) and Zang1817. In this study, we fine-mapped QFLANG-4B and validated its associated genetic effect. We developed a BC3F3 population using ND3331 as the recurrent parent through marker-assisted selection, as well as near-isogenic lines (NILs) by selfing BC3F3 plants carrying different heterozygous segments for the QFLANG-4B region. We obtained eight recombinant types for QFLANG-4B, narrowing its location down to a 5.3-Mb region. This region contained 76 predicted genes, 7 of which we considered to be likely candidate genes for QFLANG-4B. Marker and phenotypic analyses of individual plants from the secondary mapping populations and their progeny revealed that the FLANG of the ND3331 allele is significantly higher than that of the Zang1817 allele in multiple environments. These results not only provide a basis for the map-based cloning of QFLANG-4B, but also indicate that QFLANG-4B has great potential for marker-assisted selection in wheat breeding programs designed to improve plant architecture and yield.

Gonadal Transcriptome Analysis Reveals that SOX17 and CYP26A1 are Involved in Sex Differentiation in the Chinese Soft-Shelled Turtle (Pelodiscus sinensis).

The Chinese soft-shelled turtle (Pelodiscus sinensis) is an important aquaculture animal in China and exhibits growth dimorphism. Single-male cultures are often selected for higher economic efficiency. However, the mechanism of sex differentiation in P. sinensis is not well-known. In this study, a comparative transcriptome analysis of male (ZZ)- and 17β-oestradiol (E2)-induced pseudo-female (ZZ + E2)-stage embryonic gonads of P. sinensis was performed. A total of 420 differentially expressed genes (DEGs), which included 271 upregulated genes and 149 downregulated genes, were identified. These DEGs were mainly involved in several sex-related pathways, such as "ovarian steroidogenesis", "steroid hormone biosynthesis", "PPAR signalling pathway", and "metabolism of xenobiotics by cytochrome P450". In addition, 50 known and novel candidate genes involved in sex differentiation, such as the male-biased genes AMH, DMRT1, TBX1, and CYP26A1 and the female-biased genes CYP1A1, RASD1, and SOX17, were investigated and identified. For further verification, the full-length cDNAs of SOX17 and CYP26A1 were obtained. SOX17 contains a 1218-bp ORF and encodes 405 amino acids containing an HMG functional domain unique to the Sox superfamily. CYP26A1 contains a 1485-bp ORF and encodes 494 amino acids. Different expression levels of SOX17 and CYP26A1 could be detected in all the tested tissues of males and females. Notably, the expression of CYP26A1 was markedly greater in the gonads of male embryos (P < 0.05) than in those of female embryos, whereas the expression of SOX17 showed the opposite trend (P < 0.05). Taken together, the RNA-seq and qRT‒PCR results suggested potential roles for SOX17 and CYP26A1 in promoting female and male gonadal development, respectively, in P. sinensis. Our results provide new evidence for the mechanism of sex differentiation in P. sinensis.

Cytochrome GmGLY1 is Involved in the Biosynthesis of Glycitein in Soybean.

Isoflavones, the major secondary metabolites of interest due to their benefits to both human and plant health, are exclusively produced by legumes. In this study, we profiled the isoflavone content in dry seeds from 211 soybean [Glycine max (L.) Merr.] accessions grown across five environments. Broad and discernible phenotypic variations were observed among accessions, regions, and years of growth. Twenty-six single-nucleotide polymorphisms (SNPs) associated with the sum of glycitein (GLE), glycitin (GL), 6″-O-acetylglycitin (AGL), and 6″-O-malonylglycitin (MGL) contents were detected in multiple environments via a genome-wide association study (GWAS). These SNPs were located on chromosome 11 (8,148,438 bp to 8,296,956 bp, renamed qGly11-01). Glyma.11g108300 (GmGLY1), a gene that encodes a P450 family protein, was identified via sequence variation analysis, functional annotation, weighted gene coexpression network analysis (WGCNA), and expression profile analysis of candidate gene, and hairy roots transformation in soybean. Overexpression of GmGLY1 increased the glycitein content (GLC) in soybean hairy roots and transgenic seeds, while CRISPR/Cas9-generated mutants exhibited decreased GLC and increased daidzein content (DAC). Haplotype analysis revealed that GmGLY1 allelic variations significantly affect the GLC accumulation. These findings enhance our understanding of genes influencing GLC in soybean and may guide breeding for lines with high and stable GLC.

The immune cells have complex causal regulation effects on cancers

BackgroundThere was a large body of evidence linking immune cells to cancer risk. However, the causal relationship between immune cells, cancer, and what genes play an important role is unclear. MethodsIn this study, we performed comprehensive two-sample Mendelian randomization analysis (TSMR) to determine the causal relationship between immune cells and common cancers. We also performed genome-annotated multiple marker analysis (MAGMA) on immune cells causally associated with cancer to identify their relevant genes and used data summary-based MR (SMR) analysis to investigate the causal relationship between their gene expression, methylation, and cancer, and further used drug prediction and molecular docking to validate the medicinal value of the targets. Finally, reverse TSMR analysis was performed on cancer and immune cells to rule out reverse causality. ResultsAfter FDR correction (PFDR < 0.05), the results showed that 2 immune cells were associated with lung cancer risk, and 1 immune cell was significantly associated with pancreatic cancer risk. The expression of OSBPL10, CHD4, SMDT1, PHETA2, and NAGA was positively and causally related to the risk of lung cancer by SMR analysis and HEIDI test. We also found that increased expression of ANP32E decreased the risk of pancreatic cancer and that the methylation level of OSBPL10, CHD4, SULF2, CENPM, and CYP2D6 had a causal association with lung cancer. The methylation level of FCGR3A was causally associated with pancreatic cancer. The results of molecular docking indicated a strong affinity between the drugs and proteins that possessed existing structural information. ConclusionThis data-driven MR study demonstrates the causal role of immune cells in cancers. In addition, this study identifies candidate genes that may be potential anti-cancer drug targets.

NAC1 transcriptional activation of LDHA induces hepatitis B virus immune evasion leading to cirrhosis and hepatocellular carcinoma development

Our study aimed to elucidate the molecular mechanisms underlying NAC1 (nucleus accumbens associated 1) transcriptional regulation of LDHA and its role in HBV immune evasion, thus contributing to the development of cirrhosis and hepatocellular carcinoma (HCC). Utilizing public datasets, we performed differential gene expression and weighted gene co-expression network analysis (WGCNA) on HBV-induced cirrhosis/HCC data. We identified candidate genes by intersecting differentially expressed genes with co-expression modules. We validated these genes using the TCGA database, conducting survival analysis to pinpoint key genes affecting HBV-HCC prognosis. We also employed the TIMER database for immune cell infiltration data and analyzed correlations with identified key genes to uncover potential immune escape pathways. In vitro, we investigated the impact of NAC1 and LDHA on immune cell apoptosis and HBV immune evasion. In vivo, we confirmed these findings using an HBV-induced cirrhosis model. Bioinformatics analysis revealed 676 genes influenced by HBV infection, with 475 genes showing differential expression in HBV-HCC. NAC1 emerged as a key gene, potentially mediating HBV immune escape through LDHA transcriptional regulation. Experimental data demonstrated that NAC1 transcriptionally activates LDHA, promoting immune cell apoptosis and HBV immune evasion. Animal studies confirmed these findings, linking NAC1-mediated LDHA activation to cirrhosis and HCC development. NAC1, highly expressed in HBV-infected liver cells, likely drives HBV immune escape by activating LDHA expression, inhibiting CD8 + T cells, and promoting cirrhosis and HCC development.

Transcriptomic and metabolomic dissection of skeletal muscle of crossbred Chongming white goats with different meat production performance

BackgroundThe transcriptome and metabolome dissection of the skeletal muscle of high- and low- growing individuals from a crossbred population of the indigenous Chongming white goat and the Boer goat were performed to discover the potential functional differentially expressed genes (DEGs) and differential expression metabolites (DEMs).ResultsA total of 2812 DEGs were detected in 6 groups at three time stages (3,6,12 Month) in skeletal muscle using the RNA-seq method. A DEGs set containing seven muscle function related genes (TNNT1, TNNC1, TNNI1, MYBPC2, MYL2, MHY7, and CSRP3) was discovered, and their expression tended to increase as goat muscle development progressed. Seven DEGs (TNNT1, FABP3, TPM3, DES, PPP1R27, RCAN1, LMOD2) in the skeletal muscle of goats in the fast-growing and slow-growing groups was verified their expression difference by reverse transcription-quantitative polymerase chain reaction. Further, through the Liquid chromatography-mass spectrometry (LC-MS) approach, a total of 183 DEMs in various groups of the muscle samples and these DEMs such as Queuine and Keto-PGF1α, which demonstrated different abundance between the goat fast-growing group and slow-growing group. Through weighted correlation network analysis (WGCNA), the study correlated the DEGs with the DEMs and identified 4 DEGs modules associated with 18 metabolites.ConclusionThis study benefits to dissection candidate genes and regulatory networks related to goat meat production performance, and the joint analysis of transcriptomic and metabolomic data provided insights into the study of goat muscle development.

Unraveling the genetic architecture of blood unfolded p-53 among non-demented elderlies: novel candidate genes for early Alzheimer's disease

BackgroundAlzheimer's disease (AD) is a heritable neurodegenerative disease whose long asymptomatic phase makes the early diagnosis of it pivotal. Blood U-p53 has recently emerged as a superior predictive biomarker for AD in the early stages. We hypothesized that genetic variants associated with blood U-p53 could reveal novel loci and pathways involved in the early stages of AD.ResultsWe performed a blood U-p53 Genome-wide association study (GWAS) on 484 healthy and mild cognitively impaired subjects from the ADNI cohort using 612,843 Single nucleotide polymorphisms (SNPs). We performed a pathway analysis and prioritized candidate genes using an AD single-cell gene program. We fine-mapped the intergenic SNPs by leveraging a cell-type-specific enhancer-to-gene linking strategy using a brain single-cell multimodal dataset. We validated the candidate genes in an independent brain single-cell RNA-seq and the ADNI blood transcriptome datasets. The rs279686 between AASS and FEZF1 genes was the most significant SNP (p-value = 4.82 × 10–7). Suggestive pathways were related to the immune and nervous systems. Twenty-three candidate genes were prioritized at 27 suggestive loci. Fine-mapping of 5 intergenic loci yielded nine cell-specific candidate genes. Finally, 15 genes were validated in the independent single-cell RNA-seq dataset, and five were validated in the ADNI blood transcriptome dataset.ConclusionsWe underlined the importance of performing a GWAS on an early-stage biomarker of AD and leveraging functional omics datasets for pinpointing causal genes in AD. Our study prioritized nine genes (SORCS1, KIF5C, TMEFF2, TMEM63C, HLA-E, ATAT1, TUBB, ARID1B, and RUNX1) strongly implicated in the early stages of AD.

Genetic variation in two candidate genes against gastrointestinal parasites in Colombian Hair Sheep

Objective. To characterize two SNP-type genetic polymorphisms in the GLI1 (rs411868094) and IL20RA (rs419463995) genes, candidates for resistance against gastrointestinal parasites in two biotypes of Colombian hair sheep. Materials and methods. From the DNA bank of the Animal Genetics Laboratory of the University of Sucre, 167 samples, belonging to the Ethiopian (n=94) and Sudan (n=73) biotypes, were analyzed by PCR and subsequent bidirectional sequencing of two SNPs in the GLI1 (T&gt;G) and IL20RA (G&gt;A) genes. Allelic and genotypic frequencies, observed (Ho) and expected (He) heterozygosity, F index, and Hardy-Weinberg equilibrium deviations (EHW) were calculated using the GENALEX version 6.5 program. Results. For the GLI1 locus, the mean genotypic frequencies were 0.155±0.07, 0.370±0.07, and 0.475±0.07 for GG, GT, and TT, respectively. In the Ethiopian biotype, the highest frequencies of the genotype of interest (GG) were found. For the IL20RA locus, the AA and AG genotypes had similar and the highest frequency (0.465±0.03) compared to the GG genotype (0.110±0.01). The genotype of interest at this locus (AA) was the most frequent in both OPC biotypes. Conclusions. The alleles of interest associated with low FEC had a low frequency for the GLI1 gene, but a high frequency for IL20RA. The Ethiopian OPC biotype showed the highest frequencies of the genotypes of interest.

Influence of clinical and genetic factors on propofol dose requirements: A genome-wide association study.

Propofol is a widely used intravenous hypnotic. Dosing is mostly based on weight with great interindividual variation in consumption. Suggested factors affecting propofol requirements include age, gender, ethnicity, anxiety, alcohol consumption, smoking, and concomitant valproate use. Genetic factors have not been widely explored. We studied 1000 women undergoing breast cancer surgery under propofol and remifentanil anesthesia. Depth of anesthesia was monitored with State Entropy TM. Propofol requirements during surgery were recorded. DNA from blood was genotyped with a genome-wide array. A multivariable linear regression model was used to assess the relevance of clinical variables and select those to be used as covariates in a genome-wide association study (GWAS). Imputed genotype data were used to explore selected loci further. In silico functional annotation was used to explore possible consequences of the discovered genetic variants. Additionally, previously reported genetic associations from candidate gene studies were tested. BMI, smoking status, alcohol use, remifentanil dose (ln(mgkg -1min -1)), and average state entropy during surgery remained statistically significant in the multivariable model. Two loci reached genome-wide significance (P < 5×10 -8). The most significant associations were for SNPs rs997989 (30 kb from ROBO3), likely affecting expression of another nearby gene FEZ1, and rs9518419, close to NALCN (sodium leak channel); rs10512538 near KCNJ2 encoding Kir2.1 potassium channel showed suggestive association (P = 4.7×10 -7). None of these SNPs are coding variants but possibly affect the regulation of nearby genes. None of the SNPs previously reported as affecting propofol pharmacokinetics or pharmacodynamics showed association in our data. In this first GWAS exploring propofol requirements, we discovered novel genetic associations suggesting new biologically relevant pathways for propofol and general anesthesia. The roles of the gene products of ROBO3/FEZ1, NALCN and KCNJ2 in propofol anesthesia warrant further studies.

Joint transcriptomic and metabolomic analysis provides new insights into drought resistance in watermelon (Citrullus lanatus)

IntroductionWatermelon is an annual vine of the family Cucurbitaceae. Watermelon plants produce a fruit that people love and have important nutritional and economic value. With global warming and deterioration of the ecological environment, abiotic stresses, including drought, have become important factors that impact the yield and quality of watermelon plants. Previous research on watermelon drought resistance has included analyzing homologous genes based on known drought-responsive genes and pathways in other species.MethodsHowever, identifying key pathways and genes involved in watermelon drought resistance through high-throughput omics methods is particularly important. In this study, RNA-seq and metabolomic analysis were performed on watermelon plants at five time points (0 h, 1 h, 6 h, 12 h and 24 h) before and after drought stress.ResultsTranscriptomic analysis revealed 7829 differentially expressed genes (DEGs) at the five time points. The DEGs were grouped into five clusters using the k-means clustering algorithm. The functional category for each cluster was annotated based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database; different clusters were associated with different time points after stress. A total of 949 metabolites were divided into 10 categories, with lipids and lipid-like molecules accounting for the most metabolites. Differential expression analysis revealed 22 differentially regulated metabolites (DRMs) among the five time points. Through joint analysis of RNA-seq and metabolome data, the 6-h period was identified as the critical period for watermelon drought resistance, and the starch and sucrose metabolism, plant hormone signal transduction and photosynthesis pathways were identified as important regulatory pathways involved in watermelon drought resistance. In addition, 15 candidate genes associated with watermelon drought resistance were identified through joint RNA-seq and metabolome analysis combined with weighted correlation network analysis (WGCNA). Four of these genes encode transcription factors, including bHLH (Cla97C03G068160), MYB (Cla97C01G002440), HSP (Cla97C02G033390) and GRF (Cla97C02G042620), one key gene in the ABA pathway, SnRK2-4 (Cla97C10G186750), and the GP-2 gene (Cla97C05G105810), which is involved in the starch and sucrose metabolism pathway.DiscussionIn summary, our study provides a theoretical basis for elucidating the molecular mechanisms underlying drought resistance in watermelon plants and provides new genetic resources for the study of drought resistance in this crop.

Development of a mobile, high-throughput, and low-cost image-based plant growth phenotyping system.

Nondestructive plant phenotyping forms a key technique for unraveling molecular processes underlying plant development and response to the environment. While the emergence of high-throughput phenotyping facilities can further our understanding of plant development and stress responses, their high costs greatly hinder scientific progress. To democratize high-throughput plant phenotyping, we developed sets of low-cost image- and weight-based devices to monitor plant shoot growth and evapotranspiration. We paired these devices to a suite of computational pipelines for integrated and straightforward data analysis. The developed tools were validated for their suitability for large genetic screens by evaluating a cowpea (Vigna unguiculata) diversity panel for responses to drought stress. The observed natural variation was used as an input for a genome-wide association study, from which we identified nine genetic loci that might contribute to cowpea drought resilience during early vegetative development. The homologs of the candidate genes were identified in Arabidopsis (Arabidopsis thaliana) and subsequently evaluated for their involvement in drought stress by using available T-DNA insertion mutant lines. These results demonstrate the varied applicability of this low-cost phenotyping system. In the future, we foresee these setups facilitating the identification of genetic components of growth, plant architecture, and stress tolerance across a wide variety of plant species.

FOXP1 regulates the development of excitatory synaptic inputs onto striatal neurons and induces phenotypic reversal with reinstatement.

Long-range glutamatergic inputs originating from the cortex and thalamus are indispensable for striatal development, providing the foundation for motor and cognitive functions. Despite their significance, transcriptional regulation governing these inputs remains largely unknown. We investigated the role of a transcription factor encoded by a high-risk autism-associated gene, FOXP1, in sculpting glutamatergic inputs onto spiny projection neurons (SPNs) within the striatum. We find a neuron subtype-specific role of FOXP1 in strengthening and maturing glutamatergic inputs onto dopamine receptor 2-expressing SPNs (D2 SPNs). We also find that FOXP1 promotes synaptically driven excitability in these neurons. Using single-nuclei RNA sequencing, we identify candidate genes that mediate these cell-autonomous processes through postnatal FOXP1 function at the post-synapse. Last, we demonstrate that postnatal FOXP1 reinstatement rescues electrophysiological deficits, cell type-specific gene expression changes, and behavioral phenotypes. Together, this study enhances our understanding of striatal circuit development and provides proof of concept for a therapeutic approach for FOXP1 syndrome and other neurodevelopmental disorders.

Effects of nonsynonymous single nucleotide polymorphisms of the KIAA1217, SNTA1 and LTBP1 genes on the growth traits of Ujumqin sheep

Sheep body size can directly reflect the growth rates and fattening rates of sheep and is also an important index for measuring the growth performance of meat sheep. In this study, high-resolution resequencing data from four sheep breeds (Dorper sheep, Suffolk sheep, Ouessant sheep, and Shetland sheep) were analyzed. The nonsynonymous single nucleotide polymorphisms of three candidate genes (KIAA1217, SNTA1, and LTBP1) were also genotyped in 642 healthy Ujumqin sheep using MALDI-TOFMS and the genotyping results were associated with growth traits. The results showed that different genotypes of the KIAA1217 g.24429511T&amp;gt;C locus had significant effects on the chest circumferences of Ujumqin sheep. The SNTA1 g.62222626C&amp;gt;A locus had different effects on the chest depths, shoulder widths and rump widths of Ujumqin sheep. This study showed that these two sites can be used for marker-assisted selection, which will be beneficial for future precision molecular breeding.